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Abstract We present the high-resolution X-ray spectra of Mkn421 obtained in November 2003 with the RGS aboard the 
XMM-Newton satellite. This Target of Opportunity observation was triggered because the source was in a high state of activity 
in the X-ray band. These data are compared with three archival RGS observations of the same source performed in November 
and December 2002 and one in June 2003. We searched for the presence of absorption features due to warm-hot intergalactic 
medium (WHIM). We identify various spectral features, most of which are of instrumental origin. With the sensitivity provided 
by our spectra we were able to identify only two lines of astronomical origin, namely features at 23.5 A, probably due to 
interstellar neutral oxygen absorption, and at 21.6 A, which corresponds to a zero-redshift OVII Kq transition. For the latter, 
we derive an upper limit to the gas temperature, which is consistent with WHIM, of a few times 10 5 K, if the gas density has a 
value of n e > 10" 5 cm -3 . 

Key words. BL Lacertae objects: general - X-rays: galaxies -X-rays: ISM - BL Lacertae objects: individual: Mkn 421 
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The baryon density at z > 2 calculated from observations of 
hydrogen and helium absorption lines in the Lya forest (Rauch 
et al. 1997) is in good agreement with standard Big-Bang 
nucleosynthesis predictions (Buries & Tyler 1998). This is 
not true at lower redshifts: at z ~ 0, all current analyses 
indicate that, after adding the well-observed contributions, the 
local baryon density is Clb a r,obs ~ 0.007 (see e.g. Fukugita, 
Hogan & Peebles 1997), much lower than the predicted value 
n bar ,n-syn ~ 0.035 - 0.04 (Buries & Tytler 1998). 
High-resolution, large-scale hydrodynamic simulations of 
galaxy formation have been used to predict the baryon distri- 
bution at the present epoch and at moderate redshift. The main 
result of such simulations is that approximately 30%-40% 
of the baryons in the present-day universe should reside in a 
Warm-Hot Intergalactic Medium (WHIM), shock-heated to 
temperatures of 10 5 — 10 7 K. Most of these warm-hot baryons 
seem to reside in diffuse filamentary large-scale structures 
with overdensities of 10-30 and not in virialized objects such 
as galaxy groups (Cen & Ostriker 1999; Dave et al. 2001). 
Numerical simulations have demonstrated that H-like and 
He-like ions of the heavy elements composing the WHIM 
can give rise to absorption lines in the soft X-ray spectra of 
background sources. OVII, OVIII and NelX should dominate 
the relative abundance distributions of collisionally ionized 
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and photoionized+collisionally ionized gas over a broad range 
of temperatures (5 x 10 5 -10 7 K, Nicastro et al. 1999), where 
the WHIM distribution peaks (Dave et al. 2001; Fang & 
Canizares 2000). Absorption lines at 21.6 A, 18.97 A and 13.4 
A, respectively, should then be observable in the soft X-ray 
spectra of strong background sources (see e.g. Aldcroft et al. 
1994; Hellsten, Gnedin & Miralda-Escude 1998; Perna & 
Loeb 1998; Fang & Canizares 2000). 

Absorption lines such as the doublet of OVI (1032 A and 
1038 A) should be detectable also in the optical-UV band 
(Mulchaey et al. 1996; Cen et al. 2001), for temperatures 
of 1 — 5 x 10 5 K, tracing the low-temperature tail of the 
WHIM distribution. Current far-UV observations have proved 
the existence of such a low-temperature component with the 
detection of OVI absorption lines up to z ~ 0.2 (e.g. Sembach 
et al. 2000; Tripp et al. 2001). 

Soft X-ray spectra provide an even better opportunity than the 
UV spectra: the detection and the study of these components is 
needed for the proper understanding of large and small-scale 
structures in the Universe, providing independent constraints 
on cosmological parameters. 

The predicted highly ionized gas, however, has been poorly 
studied so far, because of instrumental limitations. New spec- 
trometers aboard Chandra (HRCS/LETG) and XMM-Newton 
(RGS) have increased the sensitivity and the resolution of 
the X-ray observatories slightly beyond the WHIM detection 
limit. At present, only the strongest of these systems (EW > 10 
mA) have been detected against the spectra of very bright 
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XMM-Newton RGS 


Revolution 


Obs. Id. 


Start time 


Total exposure" 


Net exposure 






Day 


Hour 


(10 4 s) 


(10 4 s) 


0532 


0136540301 


04/11/2002 


00:44:59 


2.3 


2.1 


0532 


0136540401 


04/11/2002 


07:41:43 


2.3 


2.1 


0546 


0136541001 


01/12/2002 


22:59:25 


7.1 


5.8 


0637 


0158970101 


01/06/2003 


12:48:50 


4.3 


3.0 


0720 


0150498701 


14/11/2003 


16:14:04 


5.7 


4.7 



Table 1. Log of the XMM-Newton RGS observations of Mkn 421 taken in 2002 and 2003 that were used during this analysis. 
A further observation of Mkn421 0537.0136540701 was also performed during the cooling ofRGSl which we have not used. a : 
total amount of Good Time Intervals. b : the reprocessing of these data failed; we exclude them from further analysis. 



background sources (Nicastro et al. 2002; Mathur, Weimberg 
& Chen 2002; Fang et al. 2002; Fang, Sembach & Canizares 
2003; Cagnoni et al. 2004; Nicastro 2003). Furthermore, only 
two of them were identified as signatures of the WHIM outside 
the Local Group: an absorbing system at z ~ 0.05 toward 
3C273 (Fang et al. 2002) and one at z ~ 0.01 toward Mkn 
421 (Cagnoni et al. 2002; Nicastro et al. 2003). 
High-resolution observations of Mkn 421 have already been 
performed with Chandra (Nicastro et al. 2001; Nicastro 2003) 
and with XMM-Newton (Cagnoni 2001; de Vries et al. 
2003). The [18-24] A spectrum of Mkn 421 changed during 
the two Chandra observations of 2000. In one, the source 
was very bright and no absorption or emission lines were 
detected. In the other, negative and positive residuals from 
the best-fit power-law model were observed. Nicastro et al. 
(2001) tentatively identified two different absorbing/emitting 
systems and proposed that they are intrinsic to the nuclear 
environments, becoming fully ionized - and thus transparent - 
as the source brightens. 

Further Chandra observations, performed while the source 
was in a very bright phase, allowed Nicastro (2003) to claim 
the presence of three absorbing systems located at z ~ 0, 
z ~ 0.011 and at z ~ 0.03, respectively. The first was 
identified as the WHIM inside the Local Group, while the last 
was explained as Mkn 421 intrinsic absorption. The system st 
z ~ 0.011, was interpreted as the WHIM outside the Local 
Group. 

Even the RGS spectrum of Cagnoni (2001) showed two 
absorbing systems, one inside the Local Group and one at 
z ~ 0.01, contributing an OVII Ka absorption line at ~ 21.8 
A. de Vries et al. (2003), instead, concentrated on the [22-24] 
A range of the RGS data and found evidence for a 23.5 
A interstellar neutral Oxygen (ls-2p) absorption feature. 
Furthermore, they showed a feature at 22.77 A, which they 
argued to be a non-Galactic OVI blend. 
In this paper we present the results of a XMM-Newton RGS 
observation of Mkn 421 performed as part of a Target of 
Opportunity (ToO) program to observe Blazars in high state of 
activity (see Tagliaferri et al. 2001). We triggered the pointing 
because the All Sky Monitor aboard Rossi XTE had been 
reporting high X-ray fluxes from Mkn 421 for several days. 
We compared our observations with 3 archival RGS spectra 
of Mkn 42 1 taken in November and Dicember 2002 and with 
one taken in June 2003, as part of two different calibration 



campaigns. 

In the following sections, we briefly describe the observa- 
tions and show the results of the spectral analysis performed on 
the whole RGS energy range. Then we concentrate on smaller 
energy ranges, looking carefully for the presence of absorption 
features whose reality is checked in two ways: first, by com- 
paring the Mkn 421 spectra with the smooth X-ray spectrum of 
the Crab nebula and second, by paying particular attention to 
the raw data to check for instrumental effects. Occasional hot 
or cool pixels are routine in CCD data as a result of cosmic-ray 
damage and other effects and are treated as part of normal data 
analysis procedures. Such defects are usually confined to single 
pixels and are thus significantly narrower than the instrumental 
line response that is principally caused by scattering from the 
gratings. Finally, we discuss the results proposing an identifi- 
cation for the possible non-Galactic absorption lines. 

2. The XMM-Newton observations and data 
reduction 

The XMM-Newton X-ray payload consists of three Wolter 
type-1 telescopes, equipped with 3 CCD cameras for X-ray 
imaging, moderate resolution spectroscopy and photometry 
(EPIC). Two of these telescopes (those carrying the MOS 
cameras) are also provided with high resolution Reflection 
Grating Spectrometers (RGS-1 & RGS-2), operating in the 
range [0.33-2.5] keV (5-38 A). Each RGS unit deflects half 
of its telescope beam, dispersing the striking X-ray light at a 
wavelength-dependent angle, thus providing a spectral resolu- 
tion ofE/AE ~ 100 - 500 (FWHM). After the launch, how- 
ever, failures in the read-out electronics of the CCD-7 (RGS- 
1) and CCD-4 (RGS-2), covering the [10.5-14] A and the 
[20.1-23.9] A ranges, respectively, reduced by a factor of 2 
the RGS effective area at these wavelengths. 
Mkn 421 was the target of a RGS calibration campaign in 
November and December 2002, aimed at improving the instru- 
mental performances by lowering the operating temperature. 
The benefits of the cooling manifested as a dramatic reduction 
of hot columns and flickering pixels, as well as an increase of 
the Charge Transfer Efficiency (see the movies at the XMM- 
Newton site). The RGS-1 and RGS-2 were cooled in the night 
between November 13-14 and 3-4, respectively. A third obser- 
vation of Mkn 421 was carried out on December 1, 2002. 
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A further calibration campaign was performed in June 2003. 
The observational time was fractioned in several short point- 
ings: in our analysis we included only the longest one (lasting 
- 43 ks). 

Finally, we triggered a 50 ks ToO observation on 12 November 
2003, but, because of the intense Solar activity, we had to re- 
schedule it two days later, when, according to the ASM, the 
source was still in a high state. The log of the analyzed obser- 
vations is given in Tabled We excluded from the analysis the 
data collected during the RGS-1 cooling night (14/11/2002), 
because of reprocessing failures. 

We reprocessed the data using the XMM-Newton Science 
Analysis System (SAS) 5.4.1 and the same calibration files 
used by the XMM-Newton Survey Science Centre (SSC) in 
the standard Pipeline Processing (PPS files). 
Since the XMM-Newton instruments are affected by periods of 
high background activity induced by solar flares, we extracted 
the light curves of both instruments from a background region 
of the CCD-9, which is the closest to the instrument axis and 
the most susceptible to proton events (Snowden et al. 2002). 
We then excluded the flaring time intervals (net exposures are 
reported in Table 0. After this filtering, we re-extracted the 
source and background spectra within the 95% and outside the 
98% of the PSF, respectively. Since the RGS wavelength cali- 
bration is strongly position-dependent, we fixed the source po- 
sition to the VLBI coordinates (Ma et al. 1998). Because of 
significantly better statistics, we focused our analysis on the 
first-order data only. 

3. Looking for lines in various segments of the 
spectra 

Below 0.5 keV there are still some uncertainties in the EPIC 
calibration. Thus, before looking for the WHIM signature, we 
examined the full range of the RGS spectra to derive and com- 
pare the Mkn 421 spectral shape, with the values obtained with 
EPIC (see Ravasio et al. 2004). Furthermore, we checked the 
cross-calibration between the RGS and the EPIC-PN detec- 
tors in the common energy range ([0.6-1.77] keV), where both 
instruments should be properly calibrated. These results are re- 
ported in the Appendix, where in Table IA.1I we also give the 
source fluxes. 

In this Section we shall look for the possible presence of 
faint absorption features in the RGS spectra, which could be 
the signature of the WHIM toward the source. According to the 
simulations on the WHIM chemical composition, the strongest 
absorption lines should be the OVII Ka (21.602 A in the ob- 
server frame) and the OVIII Ka (18.97 A) (Hellsten, Gnedin 
& Miralda-Escude 1998). Therefore we concentrated our anal- 
ysis on small energy ranges (2-3 A wide) centered on these 
wavelengths as well as on the [22.5-24] A range, where the 
interstellar OI ls-2p absorption line (~ 23.5 A) was already 
observed by XMM-Newton and Chandra toward Mkn 421 (de 
Vries et al. 2003) and toward other sources (e.g. PKS 2155- 
304, Cagnoni et al. 2004). 

Because of failures in the read-out electronics of the CCD^ of 
the RGS-2, covering the [20.1-23.9] A band, we focused on 
the RGS-1 data, using the RGS-2 spectra, where available, to 



check the reality of the possible features. We also compared the 
RGS-1 spectra with a RGS-1 spectrum of a powerful Galactic 
source, the Crab. 

Using XSPEC 11.2.0, we extracted the unbinned, unfolded 
spectra in the [18-20] A, [21-22.3] A and [22.5-24] A in- 
tervals and analyzed them with Sherpa 2.2. 1 which can better 
work in the wavelength space. Each feature that we studied, 10 
in total, is identified by a unique number throught the paper and 
in the figures. 

3.0. 1. Wavelength calibration 

We used the interstellar neutral oxygen feature at ~ 23.5 A, (OI 
ls-2p, ls 2 2s 2 2p 4 3 P° - ls2s 2 2p 5 3 P°) to determine the abso- 
lute line position in the RGS-1 spectra, which is fundamental 
to obtain reliable redshifts of the possible absorbing systems. 
This line was already observed toward Mkn 421 by XMM- 
Newton and by Chandra (de Vries et al. 2003) as well as toward 
other sources, such as PKS 2155-304 (Nicastro et al. 2002; 
Cagnoni et al. 2004) or H 1821+643 (Mathuret al. 2003). 
In Fig. [T] we show the [22.5-24] A RGS-1 spectra of our 
ToO observation (left) and of the archival data (right). We plot 
the spectra and the best-fit power-law models (upper panel), 
the residuals (mid panel) and the RGS-1 effective areas dur- 
ing each XMM-Newton observation (bottom panel). For the 
archival data, we plot the first and the second exposures of 
November 4, 2002 as solid and dotted lines, respectively, the 
December 1 st , 2002 data as a short-dashed line and the June 
1 st data as a long-dashed line. Residuals at ~ 23.5 A can be 
observed in the middle panels of Fig.[Q which are very likely 
produced by the interstellar neutral oxygen. Firstly, we fitted 
the ToO observation with a power-law + one Gaussian, then 
we simultaneously fitted the four archival spectra. Finally we 
fitted the five spectra together. The best-fit line position during 
the ToO observation is 23.510 ± 0.007 A, which is consistent 
with the archival data result (23.507 ± 0.013 A). Fitting simul- 
taneously all the spectra we obtained A = 23.510 ± 0.007 A. 
This value is slightly higher than the theoretical position found 
by Mc Laughlin & Kirby (1998; A = 23.467 A), but it is 
consistent with the results of many authors, obtained through 
different experimental techniques (see Table 0. In particu- 
lar, it is in good agreement (< 1 a) with other XMM-Newton 
and Chandra observations. Therefore, we shall not apply wave- 
length corrections to our data. 

3.1. The [18-20.3] A spectra 

In this energy range we can directly compare the RGS-1 and 
RGS-2 data. In Fig. we show the RGS-1 (dark gray lines) 
and the RGS-2 spectra (light gray lines) of the ToO observa- 
tion, together with the best-fit absorbed power-law model (dot- 
ted lines). We report also the residuals (in terms of sigmas, mid 
panels) and the effective area of the instruments (bottom pan- 
els) during the exposure. In Fig. [5] we show the same plots for 
the archival observations. We display as solid and dotted lines 
the first and the second exposures of November 4, 2002, re- 
spectively, as a short-dashed line the December 1 st , 2002 data 
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Figure 1. The ToO (left) and the archival (right) RGS-1 spectra of Mkn 421 in the [22.5-24] A band, together with the best-fit 
power-law models (upper panels). In the right panel, the solid and the dotted lines refer to the first and the second exposures of 
November 4, 2002, respectively, the short-dashed one represents the December 1 st , 2002 data while the long-dashed data are 
those taken in June 2003. In the middle panel we report the residuals left by the best-fit models and in the bottom panel the 
RGS-1 effective areas. 



A (A) 


Method 


Ref. 


23.467 


Theoretical 


1 


23.489 ± 0.004 


Experimental 


2 


23.536 ± 0.002 


Experimental 


3 


23.52 ± 0.03 


Chandra, X06 14+091 


4 


23.509±° : g°g 


Chandra, PKS 2155-304 


5 


23.50 ±0.01 


XMM 


6 


23.510 ±0.015 


XMM, PKS 2155-304 


7 


23.510 ± 0.007 


XMM, Mkn 421 


ToO obs. 


23.507 ± 0.013 


XMM, Mkn 421 


Arch, data 


23.510 ±0.013 


XMM, Mkn 421 


Total 



Table 2. Neutral oxygen wavelengths (ls-2p transition). 1) Mc 
Laughlin & Kirby, 1998; 2) Krause 1994; 3) Stolte et al. 1997; 
4) Paerels et al. 2001; 5) Nicastro et al. 2002; 6) de Vries et al. 
2003; 7) Cagnoni et al. 2004. 



and as a long-dashed line the June 1 st data. 
Several features with significance > 3 a are shown in all the 
RGS-1 spectra. Two of them, the numbers 1 and 4, are located 
in some of the archival observations at the same energies as 
large structures in the effective area curves. Furthermore, they 
cannot be observed in the corresponding RGS-2 residuals, in 



regions where the effective areas are smooth. Their origin is 
very likely instrumental and we shall exclude them from fur- 
ther analysis. 

We then fitted the RGS-1 spectra again, with an absorbed 
power-law model plus four Gaussian profiles to reproduce the 
residuals 2, 3, 5 and 6. In Table[5]we report the best-fit param- 
eters of each Gaussian. 

Besides investigating the corresponding RGS-2 spectra, we 
also checked the reality of these features by studying a RGS- 
1 spectrum of the Crab taken in August 8, 2002 (Obs. Id. 
0153750501). The Crab nebula, a Galactic source, has a very 
intense featureless power-law continuum and is therefore very 
useful to discriminate between the possible origins of the ob- 
served features. We reduced the Crab data as described in the 
previous Sections. In Fig. |4] we report the residuals (upper 
panel) given by the best-fit power-law model in the [18-20] 
A range and the relative effective area (lower panel). 

We summarise our results as follows: 

- Feature 2 (18.680 A): we found large residuals (5 — 6<r) 
close to the observer-frame wavelength of the OVII K/3 
transition, in a smooth region of the RGS-1 effective 
area. This identification, however, is still controversial. The 
RGS-2 data do not unambiguously confirm the RGS-1 re- 
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Figure 3. The RGS-1 [18-20.3] A (left) and the RGS-2 [18-20] A (right) archival spectra of Mkn 421, together with the best-fit 
power-law models (upper panels). The solid and dotted data refer to the first and the second exposures of November 4, 2002, 
respectively, the short-dashed ones represent the December 1 st , 2002 data while the long-dashed ones the June 2003 data. In the 
mid panel we report the residuals left by the best-fit models and in the bottom panel the effective areas. We limit the RGS-2 plot 
to the [18-20] A range because of the CCD-4 failure. 



RGS-1 [18-20] A 



Feature 


A 


FWHM 


k 


number 


(A) 


(mA) 


(cts cm -2 s _1 A -1 ) 


ToO observation 


2 


18.661 ±0.004 


24t» 


-0.0040 ±0.0011 


3 


18.937 ±0.001 


8± 1 


0.0111 ±0.0023 


5 


19.662 ± 0.005 


30 ± 10 


-0.0040 ± 0.0010 


6 


19.890 ±0.001 


8.5±i-° 


-0.0098 ± 0.0031 


Archival data 


2 


18.683 ±0.001 


15.2^ 


-0.319 ±0.048 


3 


18.9585 ± 0.001 


14.5 ± 2 


-0.349 ± 0.051 


5 


19.681 ± 0.002 


15 ±3 


-0.266 ±0.053 


6 


19.910 ±0.002 


16 ±3 


-0.303 ±0.050 


Total 


2 


18.680 ± 0.004 


33+'^ 


-0.147 ±0.026 


3 


18.953 ± 0.003 


2815 


-0.177 ±0.029 


5 


19.676 ± 0.003 


2315 


-0.166 ±0.032 


6 


19.904 ± 0.004 


33l? 


-0.128 ±0.027 



Table 3. Best-fit parameters of the power-law plus four 
Gaussians model reproducing the [18-20] A RGS-1 spectra of 
Mkn 421. The error bars refer to the 90% confidence intervals. 



suits. During the first ToO observation and during the sec- 
ond one of November 4, 2002, the RGS-2 effective area 
was smooth and we observed a very narrow feature at 
a slightly higher wavelength (18.697 ± 0.002 A, FWHM 
=14 ± 3 mA and A = 18.718 ± 0.004 A, FWHM = 17 ± 8 
mA, respectively). However, a large structure is present 
in the RGS-2 effective areas of the other observations. 
Furthermore, we also observed large residuals (~ 8 a) in 
the RGS-1 Crab spectrum (18.674+g;^ A). This implies 
a Galactic or an instrumental origin: the astronomical na- 
ture of this feature is still questionable. 

- Feature 3 (18.953 A): large residuals (5 — 6 a) are observed 
at the zero-redshift wavelength of the OVIII Ka line. This 
identification also is quite doubtful. This line is located 
close to a RGS-1 effective area feature (at ~ 18.9 A) and 
cannot be observed in the RGS-2 spectra, where the effec- 
tive areas are smooth. If the RGS-2 response is well cal- 
ibrated, this feature must be instrumental and the RGS-1 
residuals are probably caused by calibration uncertainties 
of the effective area structure at ~ 18.9 A. 

- Feature 5 (19.676 A): this ~ 3 a feature is observed in all 
the RGS-1 spectra, but it is not present in the corresponding 
RGS-2 data. Also the RGS-1 Crab spectrum displays small 
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Figure 2. The RGS-1 (dark gray data) and the RGS-2 [18- 
20.3] A spectra (light gray data) of Mkn 421 during the ToO 
observation of November 2003. We also show the respective 
best-fit power-law models (dashed lines). In the mid panel we 
report the residuals left by the best-fit models and in the bot- 
tom panel the effective areas. The RGS-2 effective area decays 
above ~ 20 A because of the CCD^I- failure. 
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Figure 4. Residuals left by the best-fit power-law model to the 
[18-20] A RGS-1 spectrum of the Crab (upper panel). In the 
lower panel we show the relative RGS-1 effective area. 



RGS-1 [21-22.2] A 


Feature 


A FWHM 


k 


number 


(A) (mA) 


(cts cm~ 2 s~ x A" 1 ) 


ToO spectrum 


7 


21.5915 ±0.013 


-0.0025 ± 0.0007 


8 


21.801 ± 0.0005 7±1 


-0.018 ±0.005 


Archival data 


7 


21.612 ±0.010 68+^ 


-0.012 ±0.002 


8 


21.819 ±0.001 12 ±2 


-0.342 ±0.068 


Total 


7 


21.603 ±0.008 81+^ 


-0.109 ±0.019 


8 


21.823 ±0.001 6±2 


-0.277 ±0.059 



Table 4. Best-fit parameters of the power-law plus two 
Gaussians model reproducing the [21-22] A Mkn 421 spectra 
taken by the RGS-1 . The errorbars refer to the 90% confidence 
intervals. 



residuals (~ 2 a) at 19.673 A. The absence of this line in 
the RGS-2 suggests a very likely instrumental origin. 

- Feature 6 (19.904 A): the RGS-1 residuals (~ 4 a) are not 
observed in the RGS-2 spectra. In this case, the comparison 
with the Crab spectrum is useless, since the corresponding 
effective area is characterized by a large feature. This line 
is probably instrumental. 

3.2. The [21-22] A spectra 

Since the ToO and the archival spectra do not display signif- 
icant features in the [20-21] A range, we now concentrate 
on the [21-22] A interval, where line detections have been 
claimed by other authors. In Fig. [5] we show the [21-22] A 
RGS-1 data of the ToO (left) and of the archival observations 
(right). Negative residuals are present at ~ 21.6 A (label 7, 
- 2 - 3 cr) and at - 21.8 A (label 8, > 2 a) both in the ToO 
and in the archival spectra. We reproduced the spectra with a 
power-law model plus 2 Gaussians. As before, we fitted the 
ToO spectrum alone, then the four archival spectra simultane- 
ously. Finally, we fitted all the spectra together. In Table |4] we 
give the two Gaussian best-fit parameters. 

Because of readout failures in the RGS-2 CCD^l we cannot 
compare the RGS-1 and the RGS-2 data. To check the reality 
of the two RGS-1 features, we investigated the spectrum of the 
Crab nebula. In Fig. [6] we show the residuals left by a power- 
law model in the RGS-1 Crab spectrum (upper panel) and the 
effective area of the instrument (lower panel). 

We summarize the results for this energy range as follows: 

- Feature 7 (21.603 A): broad residuals can be observed in 
each Mkn 421 spectrum at the observer-frame wavelength 
of the OVII Ka transition. The line is not present in the 
Crab spectrum. We believe it is caused by an astronomical 
OVII Ka absorbing system. 

- Feature 8 (21.823 A): large residuals can be observed at this 
wavelength (up to ~ 4cr, December 1 st ). Unfortunately, 
the comparison with the Crab data is useless because of a 
large effective area structure. In previous papers, this line 
was usually interpreted as an OVII Ka absorbing system 
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Figure 5. The RGS-1 spectra of Mkn 421 in the [21-22] A range. In the left panel we show the November 2003 ToO spectrum 
and in the right panel the archival data (upper panels), together with the best-fit models. On the right, we show with a solid 
and a dotted line the two exposures of November 4, 2002, with a short-dashed line the December 1 st , 2002 data and with a 
long-dashed line the June 2003 data. In the middle panels we display the residuals left by the best-fit models and in the lower 
panels the RGS-1 effective areas of the RGS-1. 




21.2 21.4 21.6 21. £ 

Wavelength (A) 

Figure 6. Residuals left by the best-fit power-law model to the 
[21-22] A RGS-1 spectrum of the Crab (upper panel). In the 
lower panel we show the relative RGS-1 effective area. 



at redshift z ~ 0.01. With the present data, however, we 
cannot draw firm conclusions about its astronomical origin. 



3.3. The [22.5-24.5] A spectra 

In Fig. G] we showed the [22.5-24.5] A RGS-1 spectra of the 
ToO (left) and of the archival observations (right). 
This energy range is characterized by the presence of broad in- 
strumental features caused by oxygen absorption as well as by 
the interstellar absorption around the oxygen K edge (see e.g. 
de Vries 2003). Combining several Mkn421 and PKS2155- 
304 RGS-1 spectra and comparing them with strongly ab- 
sorbed Galactic sources, de Vries et al. (2003) showed the pres- 
ence of instrumental features around 23.05 A and 23.35 A. 
The updated calibration files we used to analyze our Mkn 421 
data account for these structures (see the large features in the 
RGS-1 effective areas of Fig. [0, even if large residuals are 
still present at ~ 23.35 A, probably caused by uncertainties in 
calibrating the instrumental molecular oxygen absorption (de 
Vries 2003). 

We observed two features at a wavelength where the effec- 
tive areas are smooth, one of which (label 10 of Fig^ is the 
already-discussed interstellar neutral oxygen ls-2p absorption 
line at - 23.5 A (see Section 4.0.1). We fitted the [22.5-24.5] 
A spectra with a power-law + 2 Gaussians models to repro- 
duce the residuals and we report the results in Table [5] As a 
comparison, we show in Fig.0the residuals left by a power- 
law model of the Crab spectrum. In this case, the extremely 
large count rate strongly suggests calibration uncertainties and 
we therefore avoide reproducing the [22.9-23.6] A residuals. 
We summarize the results as follows: 



8 



Ravasio et al.: A XMM-RGS search for WHIM features in the X-ray spectra of Mkn421 




23 23.5 
Wavelength (A) 



Figure 7. Residuals left by the best-fit power-law model to the 
[22.5-24] A RGS-1 spectrum of the Crab (upper panel). In the 
lower panel we show the corresponding RGS-1 effective area. 
We did not display the residuals in the [22.9-23.6] A range 
which is affected by large instrumental features. 



- Feature 9 (22.778 A): each Mkn 421 spectrum display 
residuals (2 — 3 a) near a small structure of the RGS-1 
effective areas (at ~ 22.75 A). This feature cannot be ob- 
served in the Crab spectrum. In a previous paper, de Vries 
et al. (2003) tentatively identified this line as an OIV blend 
produced in the local intergalactic medium. 

- Feature 10 (23.510 A): this line was observed in all the 
RGS-1 spectra and can be identified with the interstellar 
neutral oxygen (ls-2p) absorption line (see the wavelength 
calibration Section). 

4. Discussion 

In the previous Sections we showed the presence of several ab- 
sorption features in the X-ray spectra of Mkn 421 taken with 



RGS-1 [22.5-24.5] A 


Feature 
number 


A FWHM 

(A) (mA) 


k 

(cts cm -2 s" 1 A" 1 ) 


ToO observation 


9 
10 


22.770 ± 0.001 9 ± 1 
23.510 ± 0.007 64t\l 


-0.011 ±0.003 
-0.004 ±0.001 


Archival data 


9 
10 


22.789+^ 15.5 ± 5 
23.510 ±0.013 64±ao 


-0.225 ± 0.0665 
-0.082 ± 0.029 


Total 


9 
10 


22.778 ± 0.001 6±i 
23.510 ± 0.007 71+" 


-0.304 ± 0.086 
-0.133 ±0.020 



Table 5. Best-fit parameters of the power-law plus three 
Gaussians model reproducing the [22-24] A Mkn 421 spectra 
taken by the RGS-1 . The error bars refer to the 90% confidence 
intervals. 



Feature 


A 




A/ AAfwhm 




number 


(A) 


ToO 


Archival 


Total 


2 


18.680 


780±^ 


1230+^ 


570t^ 


3 


18.953 


2370 ± 300 


1310 ± 180 


68011°° 


7 


21.603 


280±?° 


320±™ 


270 ± 50 


8 


21.823 


3110 ± 440 


1820 ± 300 


3640 ± 1200 


9 


22.778 


2530 ± 280 


1470 ± 470 


3800±f 2 ° 


10 


23.510 


370 ± 70 


3 70 +no 


330±£ 



Table 7. A/ AAfwhm ratios for the possible astronomical fea- 
tures observed in the RGS-1 spectra of Mkn 421. The instru- 
mental resolving power is A/AAfwhm < 500. 



the RGS-1 aboard XMM-Newton. The X-ray spectra are sim- 
ilar, displaying the same absorption lines during all 5 obser- 
vations over one year. We discarded some of the lines because 
they are very likely instrumental. Among the rest, we found the 
well-known interstellar neutral oxygen ls-2p line at ~ 23.5 
A (see e.g. de Vries et al. 2003). We observed features at the 
observer-frame wavelengths of the expected WHIM lines, i.e. 
at - 18.6 A (OVII K/3), at - 18.95 A (OVIII Ka) and at 
~ 21.6 A (OVII Ka), which were also reported in some pre- 
vious works (e.g. Cagnoni 2001; Nicastro 2003). We cannot 
confirm the detection of the local NelX absorption feature at 
~ 13.4 A (Cagnoni 2001; Rasmussen et al. 2003). We also 
found two small features, at ~ 21 .8 A and at ~ 22.78 A, which 
were identified in previous papers as a z ~ 0.01 OVII Ka line 
(Cagnoni 2001; Nicastro 2003) and as a local OIV blend (de 
Vries et al. 2003), respectively. 

However, the astronomical origin of some of these lines is 
doubtful. The comparison of the RGS-1 data with the cor- 
responding RGS-2 spectra (where available), with a RGS-1 
spectrum of the Crab nebula and with the literature data sug- 
gests that some of these lines are probably instrumental. The 
absence of an RGS-1 feature in the corresponding RGS-2 
spectrum strongly points towards an instrumental origin of the 
line. Also the literature data are not conclusive: none of the re- 
ported lines were always detected during the XMM-Newton 
and the Chandra observations. We resume these comparisons 
in Table|6] 

The suspicions about the astronomical nature of some of the 
lines are furtherly strengthened by their extreme narrowness 
(see Table Q. The features at 18.95 A, ~ 21.8 A and at 
— 22.78 A are much smaller (by a factor > 4) than the RGS- 
1 spectral resolution (A/ AAfwhm < 500; see XMM-Newton 
users' handbook). The feature at ~ 18.6 A is marginally con- 
sistent with the instrumental performances only during the ToO 
observation. Even in this case, however, both the corresponding 
RGS-2 line and the feature in the Crab spectrum are narrower 
than the RGS response (A/ AAfwhm > 1000), supporting an 
instrumental nature. We therefore conclude that while the fea- 
tures at ~ 21.6 A, at ~ 23.5 A can be identified as astronomi- 
cal lines, all the others are probably of instrumental origin (al- 
though some of them have been reported as real astronomical 
lines in previous works by other authors). 
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Feature 
number 


RGS-1 
A (A) 


RGS-2 


Crab 


Literature 


Other 
sources 


2 


18.680 ± 0.004 


18.697 (ToO) 
18.718(4/11/02) 
instrumental 


18.674™ 


YES (C01-XMM) 
NO (NOl-Cha) 
YES (N03-Cha) 


PKS 2155-304 (N02-Cha) 
PKS 2155-304 (C04-XMM) 


3 


18.953 ± 0.003 


NO 


instrumental 


YES (C01-XMM) 
NO (NOl-Cha) 
YES (N03-Cha) 


PKS 2155-304 (N02-Cha) 
3C 273 (F03-Cha) 
NO PKS 2155-304 (C04-XMM) 


5 


19.676 ± 0.003 


NO 


NO 


NO 


NO 


6 


19.904 ± 0.004 


instrumental 


instrumental 


NO 


PKS 2155-304 (C04-XMM) 
NO PKS 2155-304 (N02-Cha) 


7 


21.603 ±0.008 


no data 


NO 


YES (C01-XMM) 
NO (NOl-Cha) 
YES (N03-Cha) 


PKS 2155-304 (N02-Cha) 
PKS 2155-304 (C04-XMM) 
3C 273 (F03-Cha) 


8 


21.823 ±0.001 


no data 


instrumental 


YES (C01-XMM) 
YES 6 (NOl-Cha) 
YES (N03-Cha) 


PKS 2155-304" (XMM) 
NO PKS 2155-304 (C04-XMM) 


9 


22.778 ±0.001 


no data 


NO 
near 
instrumental 
feature 


YES (dV03-XMM) 


PKS 2155-304 (dV03-XMM) 
PKS 2155-304 (dV03-Cha) 
PKS 2155-304 (C04—XMM) 
NO ScoXl(dV03-XMM) 
NO 4U 0614+91 (dV03-XMM) 


10 


23.510 ±0.007 


no data 




YES (dV03-XMM) 


PKS 2155-304 (POl-Cha) 
PKS 2155-304 (N02-Cha) 
PKS 2155-304 (C04-XMM) 



Table 6. Summary of absorption features observed in the RGS spectrum of Mkn42 1 . a : archival data analyzed by us. b : detected 
in one of the two observations. C01: Cagnoni (2001); N01: Nicastro et al. (2001); P01: Paerels et al. (2001); N02: Nicastro et al. 
(2002); N03: Nicastro et al. (2003); dV03: de Vries (2003); F03: Fang et al. (2003); C04: Cagnoni et al. (2004). 



The feature located at 23.510 ± 0.007 A is the well-known 
interstellar neutral oxygen (ls-2p) absorption line which we 
discussed in a previous Section. The feature at 21.603 ± 0.008 
A is very close to the OVII Ka transition (A = 21.602 A in 
the observer frame), the strongest WHIM signature predicted 
by the simulations (e.g. Hellsten, Gnedin & Miralda-Escude 
1998). An OVII absorbing system was already observed at zero 
redshift toward Mkn 421 (Cagnoni 2001; Nicastro 2003) as 
well as toward other sources as PKS 2155-489 (Nicastro et 
al. 2002; Cagnoni et al. 2004), 3C273 (Fang et al. 2002) and 
H 1821+643 (Mathur et al. 2003), while Mckernan et al. (2004) 
studied the sightlines toward 15 AGNs, finding evidence of lo- 
cal hot gas in various cases. This has been attributed either to 
the WHIM within the local group of galaxies (e.g. Nicastro et 
al. 2002, McKernan et al. 2004) or to radiatively cooling gas 
inside our Galaxy (e.g. Heckmann et al., 2002; McKernan et 
al. 2004). 

The average equivalent width of the 21.603 A line is 
EW=9.6l3 ^ consistent with that found by Cagnoni 

(2001) in her composite spectrum (EW=12.67± mA). Since 
the equivalent width of our OVII Ka line falls in the linear 
branch of the curve of growth calculated by Nicastro et al. 

(2002) , it should be produced in an unsaturated absorption 
regime. Using Fig. 4 of Nicastro et al. (2002), we can there- 
fore calculate the column density of OVII toward Mkn 421, 
Novn ~ 4 x 10 15 cm~ 2 . A similar result is also obtained us- 
ing the curves of growth calculated by Mathur et al. (2003): the 



OVII column density is in the range [4 — 10] x 10 15 cm~ 2 , de- 
pending on the assumed velocity parameter of the gas. 
Assuming an upper limit to the equivalent width of the non- 
detected OVIII Ka, we could set an upper limit to the temper- 
ature of the absorbing gas. If a line is not saturated, as in our 
case, the equivalent width produced by an ion X % can be writ- 
ten as EW(X') oc A(X)n X i (Nicastro et al. 1999b), where 
A(X) is the relative abundance of the element X compared to 
H and n X i is the relative density of the ion i of the element 
X. We assumed therefore that the OVIII Ka at ~ 18.97 A was 
characterized by the same FWHM of the observed OVII Ka 
(FWHM = 0.081±o;^ A) and by an amplitude of three times 
the uncertainty on the continuum at the line wavelength. We ob- 
tained an average EW=1.3 ± 0.3 A. Using the OVIII/OVII ratio 
vs temperature T plot calculated by Nicastro et al. (2002), we 
obtained an upper limit for the gas temperature T < 1.6 x 10 6 
K in the case of a Galactic density (n e = 1 cm -3 ). In the case 
of the extreme extragalactic low-denisity solution, n e = 10 -6 
cm~ 3 , the gas temperature would be much lower than 10 5 K, 
while for higher density values, up to n e = 5 x 10~ 5 cm~ 3 
(Nicastro et al. 2005a), the gas temperature would be compat- 
ible with values of the order of few times 10 5 K. Heckman 
et al. (2002) suggested that the zero redshift X-ray absorption 
lines observed toward PKS 2155-304 by Nicastro et al. (2002) 
originated in a radiatively cooling gas inside our Galaxy. This 
would be the case also for Mkn 421 if we assume a low extra- 
galactic density. The gas would then be too cool to be identified 
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as the searched-for WHIM, whose temperature is predicted to 
be greter than 10 5 K (e.g Hellsten, Gnedin & Miralda-Escude 
1998). However, if we assume a larger density value, then our 
data are still compatible with a WHIM with a temperature of 
the order of a few times 10 5 K. 

5. Conclusions 

We presented the analysis of the data of the high resolution 
RGS spectrometers aboard XMM-Newton during a ToO obser- 
vation of Mkn 421 performed in November 2003. The pointing 
was triggered since the source was in a high state of activity 
in the X-ray band. We compared these data with 3 archival 
RGS observations of the same source performed in November 
and December 2002 and with one performed in June 2003. We 
summarize the main results: 

- The integrated [7-36] A RGS-1 and RGS-2 spectra are 
well fitted by a broken power-law model, with a hard spec- 
tral index (ai ~ 0.6 - 0.9) below 0.7-0.8 keV softening 
toward higher energies (a^ ~ 1.0 — 1.3). The RGS-2 spec- 
tra show systematically softer slopes and lower fluxes. 

- We compared the RGS spectra to the corresponding EPIC- 
PN spectra in the common well-calibrated energy range: 
[0.6-1.77] keV. The EPIC-PN spectra are significantly 
softer, always showing spectral indices a > 1, with fluxes 
larger by 5-10%. 

- We focused on small sections, 2-3 A wide, of the RGS- 
1 spectra, looking for absorption features. We found sev- 
eral lines which are common to all the observed spectra. 
However, we found evidence that several of them are in- 
strumental: their proximity to effective area features, their 
absence in the corresponding RGS-2 spectra or their ex- 
treme narrowness (much smaller than the instrumental res- 
olution). We found that only two features are very likely of 
astronomical origin: one at ~ 23.5 A and one at ~ 21.6 A. 

- We identified the feature at 23.5 A as interstellar neu- 
tral oxygen absorption, as proposed e.g. by de Vries et al. 
(2003). We found a best-fit wavelength of 23.510 ± 0.007 
A, slightly larger than the theoretical position proposed by 
Mc Laughlin & Kirby (1998), but fully consistent with the 
experimental results of several authors. 

- The feature at ~ 21.6 A corresponds to a zero-redshift 
OVII Ka transition. Assuming a 3a upper limit on the non- 
detected corresponding OVIII Ka, we calculated an upper 
limit for the gas temperature. We found that an extragalac- 
tic gas with a very low density would be much colder than 
the predicted WHIM temperature range. In this case this 
line would have a Galactic origin, produced by a radiatively 
cooling gas, as proposed by Heckman et al. (2002). For an 
extragalactic gas with higher density, WHIM with a tem- 
perature of a few times 10 5 K would still be compatible 
with the data. 

We conclude from our analysis that, with the current knowl- 
edge of the XMM-Newton grating performances and the sensi- 
tivity provided by our spectra, we could not find firm evidence 
of the Warm/Hot Intergalactic Medium toward Mkn 421. 
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Note added in proof After this paper was accepted, the 
Nicastro et al. (2005b) paper appeared in the literature. Nicastro 
et al. were able to observe Mkn 421 in much higher states with 
the grating on board the Chandra satellite. Thanks to the higher 
signal to noise ratio they were able to identify many more as- 
tronomical lines that allow them to derive for the first time a 
population of baryons in two intervening WHIM systems at 
z = 0.011 and z = 0.027 and to study in detail the zero- 
redshift system (Williams et al. ApJ submitted). 
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Appendix A: RGS spectral analysis and 
comparison with EPIC-PN 

First, we analyzed the RGS-1 and RGS-2 spectra in the full 
[0.34-1 .77] keV energy range ([7-36] A), where the calibration 
uncertainties are smaller than 10%. We rebinned the RGS data 
to have at least 2000 counts per bin and added a systematic 
error of 3%. We fitted the data with an absorbed power-law and 
a broken power-law model, fixing the absorption parameter to 
the Galactic value (N H = 1.61 ± 0.1 x 10 20 cm~ 2 ; Lockman 
& Savage 1995). In Table lA. ll we report the best-fit parameters 
of each observation. The RGS-2 spectra display systematically 
steeper slopes and lower fluxes than the RGS-1, as reported in 
the calibration paper on the radio-loud narrow-line Seyfert 1 
galaxy PKS 0558-508 (Kirsch 2003). 

The power-law model cannot reproduce the spectra: the \r 
are always larger than 2 (e.g. Fig. IA.lt . The broken power- 
law model significantly improves the quality of the fit: the 
spectra of Mkn 421 are hard up to ~ 0.7 — 0.8 keV and 
soften toward higher energies. We also tried to reproduce the 
steepening with the logarithmic parabolic model proposed by 
Massaro et al. (2003a,b), which we already applied to the cor- 
responding EPIC-PN spectra (Ravasio et al. 2004), as well 
as to the BeppoSAX spectra of other similar sources, such as 
1ES 1959+650 (Tagliaferri et al. 2003). This curved model re- 
produces the data well, even if, in all cases, the broken power- 
law gives better results. 

To check the quality and the good calibration of these data, 
we compared them to the corresponding EPIC-PN spectra in 
the common well-calibrated energy range ([0.6-1.77] keV). 
The RGS and the EPIC-PN data are well fitted by a broken 
power-law model softening toward higher energies. The EPIC- 
PN spectra are significantly softer, always with slopes a > 1, 
while the RGS spectra are hard up to ~ 0.8 keV. The greater 
softness of the EPIC-PN was shown also by Kirsch (2003), us- 
ing an XMM-Newton observation of PKS 0558-508. This is 
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channel energy (keV) 



Figure A.2. EPIC-PN (filled circles), RGS-1 (empty circles) 
and RGS-2 (crosses) spectra of Mkn 421 during the ToO ob- 
servation of November 14, 2003, fitted together with a broken 
power-law model (with Galactic absorption value). For plot- 
ting purposes, we rebin the the RGS data to have at least 15000 
counts in each bin. The RGS spectra are significantly harder 
than the EPIC-PN spectrum. For clarity, we omitted the resid- 
ual error bars. 



very clear in Fig. IA.2l where we plot the RGS and the EPIC- 
PN spectra of the ToO observation. Furthermore, the EPIC-PN 
fluxes are larger by ~ 5 — 10%. As an example, we report in 
Table IA.2I the best-fit parameters of the power-law and bro- 
ken power-law models for the November 2003 ToO observa- 
tion. These calibration differences, however, are not important 
for our search for WHIM features. We focused on such small 
sections of the RGS spectra (2-3 A wide) that the continuum 
uncertainties can be neglected. 
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Figure A.l. The RGS-1 spectrum of the ToO observation of November 14, 2003. An absorbed power-law model is not able to 
reproduce the data well (left panel), while a broken power-law model fits the data better (right panel). We fixed the absorption 
parameter to the Galactic value. The spectra were rebinned in order to have at least 2000 counts per bin. The other RGS spectra 
display similar behaviors. 
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-2 
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-1 
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3.16 
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0.691^4 


1 20+ 003 
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3.19 
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-2 


0.94 
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3.05 
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-1 
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-1 


0.64 ±0.02 
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1.92 
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-1 
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2.3/139 


RGS- 


-2 
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1.86 
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-1 
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-1 
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4.87 
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-2 
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4.52 
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Table A.l. Best-fit parameters of the RGS spectra reproduced with power-law, broken power-law and parabolic models. We 
fixed the absorption parameter to the Galactic value. The error bars refer to the 90% confidence interval. a : erg cm~ 2 s _1 . 
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Inst. 
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(keV) 
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RGS- 
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6.6/45 
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1.02 ±0.06 
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Table A.2. Best-fit parameters of the ToO RGS and EPIC-PN spectra fitted with power-law and broken power-law models, in 
the common energy range ([0.6-1.77] keV). The absorption parameter was fixed to the Galactic value. Systematic errors of 3% 
and of 0.5 % were added to the RGS and to the EPIC-PN data, respectively. The error bars refer to the 90% confidence interval 
for one parameter. a : erg cm~ 2 s _1 . 



